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THEORETICAL STOKY OF THE TRANSONIC LIFT OF A DOUBIE-WEIGE 
PROFILE WITH DETACHED BOW WAVE 1 


By Walter G. "Vincent! and Cleo B. Wagoner- 


SUMMARY 


A theoretical study is described of the aerodynamic characteristics 
at small angle of attach of a thin, double -wedge profile in the range 
of supersonic flight speed in which the bow wave is detached. The 
analysis is carried out within the framework of the transonic (nonlinear) 
small -disturbance theory, and the effects of angle of attack are regarded 
as a small perturbation on the flow previously calculated at zero angle. 
The mixed flow about the front half of the profile is calculated by 
relaxation solution of a suitably defined boundary-value problem for the 
transonic small -disturbance equation in the hodograph plane (i.e., the 
Tricomi equation) . The purely supersonic flow about the rear half is 
found by an extension of the' usual numerical method of characteristics. 0 / 
Analytical results are also obtained, within the framework of the same 
theory, for the range of speed in which the bow wave is attached and the 
flow is completely supersonic. 

The calculations provide, for vanishingly small angle of attack, 
the following information as a function of the transonic similarity 
parameter: (1) chordwise lift distribution, (2) lift-curve slope, and 
(3) position of center of lift. As in previous studies, the aerodynamic 
characteristics of a profile of given thickness ratio show little varia- 
tion with free-stream Mach number as the Mach number passes through 1. 

As the Mach number is increased to higher values, however, the lift- 
curve slope rises to a pronounced maximum in the vicinity of shock 
attachment and then declines. Correspondingly, the center of lift moves 
forward toward the leading edge and then returns aft. These findings 
are in marked contrast to the behavior of the drag coefficient at zero 
angle of attack, which was found in earlier work to decrease monotoni- 
cally as the Mach number increased above 1. At Mach numbers above that 
for shock attachment, the results of the present calculations tend toward 
those given by classical linear theory. 

"’"Portions of this work were reported at the VTIIth International 
Congress on Theoretical and Applied Mechanics, Istanbul, Turkey, 

August 20-28, 1952. 
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INTRODUCTION 


The theoretical problem of the transonic flow over a thin, double- 
wedge profile at zero angle of attack has been treated in several papers 
in recent years. These papers have in common that they employ the 
simplifying concepts of the transonic smaU -disturbance theory and 
utilize the holograph transformation to linearize the resulting mathe- 
matical problem. Following this approach, Guderley and Yoshihara 
(reference l) began by solving the problem for a free -stream Mach 
number of 1, using analytical methods for the mixed flow over the front 
wedge and the method of characteristics for the purely supersonic flow 
over the rear. Somewhat later, the present authors, using a combina- 
tion of relaxation methods and the method of characteristics (refer- 
ences 2 and 3) , extended the results to free-stream Mach numbers greater 
than 1, where a detached bow wave occurs ahead of the profile. At about 
the same time. Cole (reference obtained an analytical solution for 
the flow over the front wedge at subsonic flight speeds, utilizing, in 
effect, the special assumption of a vertical sonic line from the 
shoulder of the wedge. More recently. Trilling (reference 5) has been 
able to remove this special assumption and, with the aid of less strin- 
gent approximations regarding the flow over the rear wedge, to extend 
the solution for the subsonic case to include the complete profile. As 
a result of these investigations, the problem of the double-wedge 
profile at zero angle of attack may be regarded as substantially solved 
within the limitations of the transonic small-disturbance theory. The 
experimental studies of Liepmann and Bryson (references 6 and 7) and 
Griffith (reference 8) indicate that the theoretical findings are in 
fundamental agreement with the physical facts. 

In a recent paper (reference 9 )? Guderley and Yoshihara have con- 
tinued their investigations of the double -wedge profile at Mach number 1 
by considering the influence of a vanishingly small angle of attack. 

The basic idea in this later work is to regard the effects of angle of 
attack as a first-order perturbation on the nonlinear flow previously 
calculated at zero angle. This approach leads to a linear boundary- 
value problem in both the physical and hodograph planes. The calcula- 
tion for the front wedge is still carried out, however, in the hodograph 
plane, since the basic procedures can then be taken over directly from 
the previous work. By this means, Guderley and Yoshihara obtain results 
for the lift-curve slope of the profile at zero angle of attack and for 
the corresponding distribution of lift along the chord. 

The aim of the present paper is to extend the results for the 
double wedge at angle of attack to the case of supersonic flight with 
detached bow wave. The fundamental ideas of Guderley and Yoshihara are 
followed in reducing the calculations for the front wedge to a perturba- 
tion problem in the hodograph plane. The detailed formulation of the 
problem is, however, necessarily different in the present case. The 
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■boundary conditions for the problem appear in terms of the results 
already obtained at zero lift (references 2 and 3 ), and the solution is 
carried out by numerical methods which differ only slightly from those 
devised for the earlier work. The lift on the rear wedge is calculated 
by an extension of the method of characteristics. The body of the 
paper is devoted to the detailed formulation of the boundary- value 
problem in the holograph plane and to a discussion of the final results. 
Noteworthy differences between the numerical procedures used in the 
present work and those already described in reference 3 are treated in 
appendices at the end of the report. 


NOTATION 
Primary Symbols 


a* critical speed (i.e., speed at which the speed of flow and 

the speed of sound are equal) 

b numerical coefficient 

(See equations (39) and*(4o).) 

c airfoil chord 


c l 


lift coefficient 


( lift per unit span ) 

\ J 


~m 


c p 


moment coefficient for moment 
( moment per unit span 


edge 


V 


% c ‘ 


) 


taken about leading 


pressure coefficient 



I w integral defined by equation (45) 

k 2 numerical constant 

(See equation (10) . ) 

M Mach number 

m slope of segment of Mach line in characteristics net 


P 


static pressure 
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Ap 

1 


local lifting pressure (i.e., difference "between static 
pressures on upper and lower surfaces) 

dynamic pressure 


t 


airfoil thickness 


V speed of flow 

x,y Cartesian coordinates 



generalized Cartesian coordinates 
(See equations (43).) 

chordwise position of center of lift 


slope of curve of lift coefficient versus true angle of 
attack evaluated at zero angle 


slope of curve of moment coefficient versus true angle of 
attack evaluated at zero* angle 


a 

P 

7 

A 


normalized angle of attack; also denotes true angle of 

dc^ 

attack when used in derivatives such as , etc. 

da 

absolute value of t| at left-hand limit of lattice 
ratio of specific heats (1.4 for air) 
basic lattice interval 


A# 


function of t] and 0 
(See equation (A 6 ).) 

normalized speed of flow 

(See equation (la) and page 10.) 




0,0 


®w 


special values of t] 

. (See sketch (m) on page 46-) 

normalized inclination of flow; 0 also denotes true inclina- 
tion of flow in equation (lb) 

(See equation (lb) and page 10.) 

normalized half angle of wedge 


^o 


transonic similarity parameter 
(See equation (13).) 
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P 

* 

6 *'a 

a,t,c 

A,B 

o 

s 

etc. 

* 

f) 

( )’ 


fluid density- 
stream function 

incremental values of stream function 
(See equations (A9) and (All).) 


Subscripts 


points in characteristics net 
(See page page 43*) 

components of total stream function 
(See equation (39)0 

conditions in free stream 

singular solution 

(See equations (a 6) and (AT)0 

value at a prescribed lattice point 


conditions at critical speed 

Superscripts 

\ 

quantities determined at zero angle of attack 

derivative with respect to normalized Angle of attack 
evaluated at zero angle 
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BOUKDARY-VALUE PROBLEM IK HODOGRAPH PLANE 


Description of Plow Field 


Sketch (a) is a drawing of the idealized, inviseid flow -which 


Shock wave 
Streamline 
Sonic tine 


M 0 >i 



Expansion \ Mach 
Compression j lines 


Sketch (a) 

may he expected about a wedge profile -when the angle of attack is 
stiff iciently less than the semiapex angle of the wedge. Sketch (b) 
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shows the corresponding hodograph representation of the flow over the 
front wedge, which is the region of prime theoretical concern. Except 
for the substitution of the detached how wave in place of the infinite 
free stream, these representations follow the lines assumed by Guderley 
and Yoshihara in reference 9* The corresponding drawings for zero angle 
of attack, which are fundamental to the present case, have been described 
in detail in references 2 and 3- 

In the present example, the path of the central streamline in the 
physical and hodograph planes is briefly as follows: The streamline 

leaves the bow wave in the physical plane (or the shock polar in the 
hodograph plane) at point A. It then proceeds with decreasing subsonic 
speed to a stagnation point 0 on the underside of the profile. At 0 
the streamline branches. The lower branch runs downstream along the 
lover surface of the profile with fixed inclination but increasing 
speed. The sonic speed is reached at the shoulder L, where the speed 
then increases diseontinuously in accord with the Prandtl -Meyer rela- 
tions. The shoulder itself maps in the hodograph onto the upgoing 
characteristic LM. The upper branch of the central streamline proceeds 
from 0 upstream along the surface of the profile. The inclination here 
is again fixed by that of the surface, and the speed increases to the 
sonic value at the leading edge J. At this point the flow is charac- 
terized by another Prandtl -Meyer expansion to supersonic speed. 

The flow configuration which should be assumed on the upper surface 
near the leading edge is open to conjecture. Since the geometrically 
available angle of turn will, for any thin airfoil, be greater than the 
130 ° permissible for expansion to a vacuum, a region of separation is 
to be expected. If the angle of attack is not too great, this region 
will probably be closed, with the central Btreamline reattaching to the 
upper surface a small distance behind the leading edge. This reattach- 
ment will be followed by a compression of the flow through a system of 
shock waves whose arrangement' is sketched only formally in the physical 
plane (and not at all in the hodograph plane, where the correct repre- 
sentation would probably lie on several sheets) . The effects of the 
flow near the leading edge will be mentioned later, but the exact 
process will remain undefined. Whatever the details, the speed on the 
upper surface will return to a subsonic value at some point K just 
downstream of a terminating, normal shock wave. From. K the central 
streamline continues at fixed inclination downstream along the upper 
surface, the speed increasing once more to the sonic value at the 
shoulder B. At this point another expansion takes place, similar to 
that which occurs at the corresponding' point on' the lower surface. In 
this case the shoulder is represented in the hodograph by the down- 
going characteristic BG. ■ 

The supersonic expansion fan from the shoulder at B (and similarly 
at L) is discussed in detail in references 2 and 3- Suffice it here to 
say that the supersonic flow field, of which the expansion fan is the 
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initial part, is separated into two regions "by the Mach line GE, which 
runs from the shoulder to the sonic point on the' how wave. (This line 
was temled the "separating" Mach line in reference 2.) The supersonic 
flow in the region upstream of the Mach line GE is interdependent with 
the subsonic field between the bow wave and sonic line. To obtain a 
solution for the front wedge, a problem in transonic flow must there- 
fore be solved for the subsonic field and the interdependent portion of 
the supersonic expansion fan. Conditions in the supersonic flow down- 
stream of the Mach line GE have no Influence upon the subsonic field. 
The continuation of the flow beyond GE can be accomplished by purely 
supersonic methods once the solution of the transonic problem is known. 

Aside from the obvious lack of symmetry in the present case, the 
main difference between the flow here and that previously studied at 
zero angle of attack is the existence in the present problem of the * 
localized supersonic region in the vicinity of the nose. As pointed 
out, conditions in this region are difficult to formulate. The problem 
has been considered by Guderley and Yoshihara (reference 9) in "the 
course of their work at Mach number 1. They find that, if the nose 
region is disregarded in the hodograph and the boundary condition 
along KB is fulfilled all the way in to 0, then the influence on the 
lift of the resulting fictitious flow at the nose is of somewhat higher 
than the second order in the angle of attack. This suggests that the 
effects of the real flow at the nose may be neglected in a first-order 
analysis such as the present. In the work' which follows, as in the 
calculations of Guderley and Yoshihara, the supersonic region at the 
leading edge will therefore be disregarded. 


Formulation of Boundary-Value Problem 


As in reference 3, the analysis is based on the equations of the 
transonic small-disturbance theory with the stream function + as the * 
dependent variable. If the effects of the flow at the nose are ignored, 
the problem of the wedge at angle of attack a is then readily formu- 
lated as a boundary-value problem in the hodograph plane. To solve this 
problem for vanishingly small d, it will be assumed that the solution^ ^ 
at angle of attack can be expressed as the sum of the basic solution 
previously obtained at zero angle plus a perturbation term aty* , where V 
is a function which does not itself 1 involve a. By consideration of the 
difference between the boundary-value problems for and $, a problem 
for the perturbation function ty’ can be formulated. The boundaries for 
this problem turn out to be the same as those for the problem at zero 
angle, and the boundary values themselves appear in terms of ^ . The 
details of these matters will now be given. The reader who is interested 
only in the results can turn directly to the section Chordwise Distribu- 
tion of Lift on page 28. 
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Basic equations . - The basic equations will be taken in the form 
given in reference 3 , that is, in terms of small disturbances from the 
critical speed a*. 2 The independent variables are the normalized speed T] 
and the normalized inclination 0 as defined by the relations 


^ Y/a *~ 1 
V 0 /a*-l 


(la) 


9 



(V 0 /a* - 1) 3/2 


(lb) 


where 

Y local speed of flow 

9 local inclination of flow relative to direction of free stream 
V 0 free -stream speed 

a# critical speed ( i . e . , speed at which the speed of flow and the 
speed of sound are equal) 

7 ratio of specific heats 

Use of these variables is equivalent to introducing the rules for tran- 
sonic similarity. In terms of the foregoing hodograph variables, the 
differential equation for the stream function y as given by the tran- 
sonic small-disturbance theory is 

" 2? i *§? = 0 (2) 


2 As discussed in several recent papers (e.’g., references 10 and ll), the 
theory can also be formulated in terms of disturbances from the free- 
stream speed V Q . This latter, less restrictive formulation reveals 
clearly the relationship which exists between the transonic small- 
disturbance theory and the familiar linear theory of subsonic or super- 
sonic flow. As shown by Spreiter, (see page 9 of reference ll), an a* 
analysis will yield values of the pressure coefficient identical to 
those of a V Q analysis provided the pressure coefficient and similar- 
ity parameter in the former case are taken as in equations (^) and (13) 
below. If this procedure is followed, the results of the a* analysis 
may even be expected to tend toward those of linear theory as the free- 
stream Mach number increases or decreases from 1. (An analytical 
example of just this behavior has been given by Bryson in appendix A 
of reference 7.) It appears, therefore, that the a^ formulation, 
when suitably used, gives results of wider theoretical validity than 
would be anticipated on the basis of its own rather restrictive under- 
lying as sump' t ion. 
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This is essentially the linear differential equation first studied "by 
Tricomi (reference 12). It is elliptic for rf<0 (subsonic speeds) 
and hyperbolic for ?f>0 (supersonic speeds). 

The transformation from the hodograph to the physical plane is 
governed by the differential relations 


= — 5 — J + v 10 


ay = _L_ (i^atf + *~d?) = d* (3b) 

m 3 P* a* 


■where x=x(rf,§ r ) and y=y(?f,<^) are physical coordinates (horizontal and 
vertical, respectively) corresponding to a given velocity 9\,6. The 
symbol p # denotes the fluid density at the critical speed a*. Within 
the approximation of the transonic small -disturbance theory, the pressure 
coefficient Cp = (p-p 0 )/q 0 can be calculated from the relation 


Cp » -2 = -2(V 0 /a*-l)(q-l) 

a * 




The local Mach number is related to the speed of flow by the equation 


M g -1 _ V_ 
7+1 a * 


(5) 


For simplicity of notation, the tilde "will be omitted from the 
symbols rf and 9 in the remainder of the development. It is to be 
understood, unless stated otherwise, that the quantities t) and 9 are 
themselves the normalized quantities defined by equations (l) . 

Problem at zero angle of attack . - When the angle of attack is zero, 
the localized region of supersonic flow at the leading edge disappears 
from sketch (a), and the flow field becomes symmetrical about the chord 
line. The corresponding boundary-value problem in the t),0 plane has 
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"been set forth in reference 3- It is restated in sketch (c), ■where hoth 
the upper and lower halves of the flow field are now included. In this 
representation, the surfaces of the wedge appear as the semi-infinite 
horizontal lines OB and OL, and the subsonic portion of the shock polar 
appears as the curve HAE. 


0 



Sketch (c) 

If the stream function for zero angle of attack is denoted hy 
1 ir=4r( T ],0) , the differential equation to "be satisfied here is given hy 
equation (2) as 

♦-in ‘ 2r <*ee “ 0 (6) 

The requirement that the flow shall he tangent to the surfaces of the 
wedge provides the boundary conditions 


^(r|,±0 w ) = 0 for T}<0 


(7) 
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where 0 W denotes the normalized half -angle at the leading edge. The 
stagnation point at the leading edge is represented in the present 
theory by the condition that 




for T ] — > -co, -Q v < Q < 0 W 


(8) 


Along the shock polar ME, the relations for an oblique shock wave 
require that 


* 


T) 


_ 1+7TJ 
+ 3+5-n 


Vl+TJ = 0 


for 


L 


(9) 


e = ± (i-tj) a/I+tJ , -i<n<o 


Along the sonic line, boundary conditions are prescribed which represent 
the influence exerted on the subsonic field by the interdependent portion 
of the supersonic expansion fans. On the basis of the procedures given 
in reference 3, this influence can be represented completely by the 
requirement that 



[±{e x -e)f /ai 


de x 


= 0 


(10) 


■where the upper signs apply for 1<0 <©v an ^- the lower signs 
for -9\r<9 -1- The constant k 2 which appears here is given by 


2 4 / 3 * 

3 a/6 r s (i/3) 


= 0. 3429 


where r(l/3) is the gamma function of the argument l/3- The use of 
the relations (10) as boundary conditions along the sonic line reduces 
the transonic problem of the flow over the front wedge to a purely 
elliptic problem in -the holograph plane. 


In addition to the foregoing conditions, a further condition is 
necessary to assure that the solution for will give the proper scale 
when transformed to the physical plane. This is furnished, for example, 
by the following expression for the half -chord of the profile, found by 
integrating equation (3a) over either OB or OL: 
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c = _2 

2 P* a * 


(7+1) (7 0 /a»-l) 
2 




(Ti,i0 w ) dt| 


(U-) 


If the chord of the profile is given, this condition, together with the 
previous conditions (7) through (10), is sufficient to determine a unique 
solution to the problem.. 

It is obvious from the nature of the boundary -value problem (and 
also from considerations of symmetry in the physical plana) that the 
solution for ijr must 'be antisymmetric with respect to 9. The problem 
can be simplified, therefore, by discarding the lower half of the holo- 
graph and replacing it by the condition 


f (tj,0) = 0 for ti <£— 1 


( 12 ) 


The resulting problem is readily solved with numerical methods by 
assuming an arbitrary value of ^ at some point (as, for example, 
point E), solving for t in the upper half of the holograph subject to 
ffche conditions (7), (8), (9), (10), and (12), and then adjusting the 
(solution to satisfy condition (ll). 


j It is apparent from the boundary conditions that the solution of 
the foregoing problem will depend on the value of the parameter 0 W , 
which defines the position of the upper and lower boundaries in the 
j holograph . This parameter is related to the more familiar transonic 
j similarity parameter ^ by the relation 


Mq 5 -1 = 2 1 ^ 3 

[(y+DCt/o )] 2/a e v s/a 


(13) 


where t/c is the thickness ratio of the complete double-wedge profile. 
In references 2 and 3 the solution of .the foregoing problem has been 
carried out for four values of 9y. 
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Problem at angle of attack. - If the supersonic region at the 
leading edge is ignored, the boundary- value problem for the wedge at 
angle of attack appears in the t |,0 plane as shown in sketch (d). 


• 8 



Sketch (d) 

The primary difference between this and the previous sketch is that the 
lines OB and OL, which represent the surfaces of the wedge, have each 
been displaced downward by an amount a, where a 1 is the angle of 
attack normalized in the same manner as the other angles of inclination 
(cf. equation (lb )). 3 ______ 

3 In reference 9 , Guderley and Yoshihara find it convenient to obtain the 
angle of attack by holding the profile fixed and changing the inclina- 
tion of the free stream. This procedure, if applied in the present 
case, would require the eventual calculation of the second derivatives 
of if on the shock polar. The present procedure, which holds the free 
stream fixed and changes the attitude of the profile, requires the cal- 
culation of only a first derivative of if at the surface of the wedge. 
Since the accuracy of numerical differentiation decreases with increas- 
ing order of the derivative, the present approach is to be preferred in 
a numerical analysis. _ 
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The stream function at angle of attack will he denoted here 
"by \Jj=i}r(Ti,0;a), the latter notation being used to indicate the depend- 
ence of t upon the parameter a. The function ^ must satisfy the 
differential equation (2), which is now written 


^T]t) - 2T i^ee 0 


(ib) 


The boundary conditions at the surface of the wedge now require that 


t(il,±0v“ a J a ) = 0 for *1^0 


(15) 


while the condition at the leading edge becomes 

>0 for > -cCj -0 v -a< g < g v ~a. 


(16) 


The shock polar ME is unaltered from the previous problem, and the 
condition on this boundary has the same form as before. The conditions 
along the segments BE and LN of the sonic line are now 


* n (O,0;a) + k. 


±0 w - a 


teCOj&jja) 

[±(e x -e )] 2/3 


d0i = 0 


(IT) 


where the lower limit of the integral has been changed in accord with 
the displacement of the points B and L. The upper signs in equation (17) 
now apply for 1 ^ 9 ^ 0 y-a and the lower signs for -0 w -a<^0 - 1. An 

expression for the half -chord of the profile can be found again by inte- 
grating equation (3a) over the line OB or OL, which gives 


£ 

2 


P* a *L 


(7+1) (V /a,, -1) 1 1/2 P° 

f J dT1 


(18) 


If the chord of the airfoil is specified - say the same as at zero angle 
of attack - then the foregoing conditions . are sufficient to determine a 
solution. No simplification based on symmetry considerations is possible 
in the present case. 

Perturbation problem . - The problem of the preceding section conceiv- 
ably could be solved by numerical methods - though with great labor - for 
arbitrary values of a. Efforts in this direction would hardly be justi- 
fied, however, in view of the fundamental omission of the localized 
supersonic flow at the leading edge. It is more reasonable to examine 
the problem for vanishingly small a, where this omission is valid and 
where there is hope that the amount of labor might be reduced. 
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To proceed along these lines, it is assumed that ty(t],0;a) may he 
expanded in a power series of the form 


t(T],G;a) = *( q,0jO) + a t a (n,0;O) + OCa 2 ) 


where, for present purposes, only terms to order a need he retained. 

The first term on the right represents the solution at a = 0 and is 
thus identical with the function' ?(t],0) previously introduced. The 
second term will he abbreviated by means of the notation ilr*(T],0)=iJr a (Ti,0;O). 
If terms of C^a, 2 ) are discarded, the expression for ^ can then 
be written 


Tlr(ib0;a) = + a (b,®) (!9) 

By comparison of the previous boundary-value problems for \|r and a 
problem for the perturbation function will now be formulated. 

The dif ferential equation for f follows at once from the differ- 
ential equations (6) and (1*0 and the substitution ( 19 ) • It is obviously 
of the same form as the previous equations, that Is, 


*’ W - 2 1 *’99- 0 (20) 

The boundary conditions appropriate to the surface of the wedge are 
established as follows: The boundary condition (15) for is first 

rewritten, with the aid of the substitution ( 19 ) , in the form 


^( T l9 ±0 W a ) + a (^,±0^) = 0 (21) 

By expanding in Taylor’s series about the lines 0 = ±0-^, the functions 
\jr and ljr’ can be written 

iKT},±0w-a) = t(q,±0w) - a ^(ti,*^) + 0(cc 2 ) (22a) 

t , (’l, ± 0w“ a ) = ’K'H**®*) - a- t’ 0 (T),±0 w ) + 0(a 2 ) (22b) 

If these 'expansions are substituted into equation (21) and ty(T),*0 w ) set 
equal to zero in accord with the boundary condition (7), one obtains 
finally for van! shingly small a 

VC 1 !**®*) = tgC 1 !**®*) tor ’iS 0 


(23) 
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This is the boundary condition for ty’ appropriate to the surface ^ of 
the profile. It -will he noted that the condition is applied in the 
holograph at the original, undeflected location of the surface 
("i.e., 0=±0 W ) • The condition depends for its application on a knowledge 
of the basic solution ty. 

The boundary condition for t' at the leading edge follows 
directly from the conditions ( 8 ) and ( 16 ) . It is the same as the corre- 
sponding condition for t, that is, 

t * — >0 for t) — >-», -dy^e^e v (2^) 

\ 

As was indicated, the functions i{r $ both satisfy the same 

linear, homogeneous boundary condition on the shock polar. It follows, 
as in the case of the differential equation, that the condition for ijr 1 
on the polar is again the same, that is, 

t' + ^^ltn t|r« = o ' (25) 

m 3+5ti y 


for 

0 = i(l— T|) a/i+t], -l^Tl^O 


The treatment of the boundary condition along the sonic line is 
complicated by the fact that the parameter a appears in the condi- 
tion (17) as a term in the lower limit of the integral. For simplicity, 
the details will be confined here to the upper segment BE of the sonic 
line. For this segment, condition (17) becomes, after substitution 
from equation ( 19 ), 


yo,e) + a.V^0,d) + *zj 


6 %( o , e x ) 


9y-q. 


(01-9) 


2/3 


dS-j_ + 




e ^ (0,6-0 

Z/3 ^ = 0 




(0i-e) 


( 26 ) 


applicable for 1^0^0 w -a. To simplify this equation, the first 
integral is rewritten 



0 w -a 


d0 = r e *e&i fiL 
(0i-0) 2/3 d01 'L (01 -0) 2/3 


■'w 


,6y-a ti0,ex) 

— - — d 0 , 

(, ex - ef /a 


( 27 ) 
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It can "be shown from Guderley’ s analysis of flow at a convex corner 
(reference 13) that, for vanishingly small values of (9 V ~&) , the 
variation of if along the sonic line must he of the form 

f(o,e) = c (e w -e) 4/s (28) 


where C is a constant for any given value of 0 W . Differentiating 
this relation, one obtains 

f e (O,0 1 ) ~ (©w-Sr) 173 


Substitution of this result into the second integral of equation (27) 
yields the fact that this integral must be proportional to a 4 / 3 . The 
first integral in equation ( 26 ) can thus be written 


f 


9 if e (O,0 1 ) 


0 w -a 


(01-0) 


2/3 


n9 * fl (O,0,) 

“i * / "-- - 57 5 + 0 <“ 4/3 > 

•' (Si-e) a/3 


(29) 


“w 


The second integral of equation (26) can be treated similarly by first 
rewriting it as 


f 

9 \ t ~ a 


*' e (oA) 

(e.-ef' 3 


d^ 


re VetOA) ^ 

V «>i-ef /3 1 ^ (%-ef /3 


dflj. 


(30) 


To deduce the variation of *' for vanishingly small (0 w -0) , it is 
first noted that a result similar to equation (28) must also hold for 
the variation of * relative to the displaced location of the shoulder, 
that is, 

*(0,0) = C (©v-a-0) 473 


The quantity C = C(a) is a differentiable function of a which reduces 
to C when a = 0. Since a will eventually be made less than any 
assignable value of (0-^-0), this expression may be expanded in the form 


*(0,0) = C (0 w -0f /3 



a 

0w-0 


+ 0(a 2 ) 


(31) 
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Now it follows from the definition of ty* that 

, ( 0 e) = lim »(0,9) - *(0,6 ) 
Y v * ' a — >0 


Substitution from equations (28) and (31) thus gives for the variation 
of ijr* in the vicinity of the shoulder 


+»(o,e) = 


lim 
a — >0 


(e, r e ) 4/3 _ ^ c (e w -e) 1/3 + o(a) 


a 


or 


t*(o,0) = c» (0 v -0) 4/3 - | c (0w-e ) 1/3 


( 32 ) 


where C’ = C a (o). This means that for vanishingly small (0 w -0) 

r Q ( 0,9 x ) - {e^-e x y z/a 

On the basis of this result equation (30) can he written 


r 

9 w -a 


r e ( 0 ,e x ) 

__ d0 

(9^9) /S 


p e r e (o,9 x ) 




+ O(aVS) 


(33) . 


If equations (29) and. (33) are substituted into equation (26) and the 
boundary condition (10) is taken into account, one then obtains for 
vanishingly small a 

r^t'aCo.ei) 

♦V ° ’*)+** J m 

0 W 

where 1^0^0 W . The boundary condition for 'i|f* along the upper segment 
of the sonic line is thus the same as the condition for The same 

result can be shown to hold along the lower segment. 


It remains to impose the condition that the chord of the airfoil 
must remain unaltered during change in angle of attack. To express this 
condition in terms of equation ( 19 ) is first substituted into 
equation (l8) to obtain 


c 

2 


2 

P* a * 


(7+1) (V 0 /a*-l) 
2 



t e ( Tj^v-a) 


+ a 


(35) 
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As in the treatment of the "boundary conditions along the upper and. 
lower boundaries, Taylor’s expansion gives 

= *e(Tb±Sw)- ct ^ 00 (ti,±6w) + OCa 2 ) . (36a) 

V e (n,*e^<i) = V e (T\,±ew)-*y e e(T\,±Ov) + o(a 2 ) ( 36 h) 

It canjbe inferred, directly from the boundary-value problem for ijr 
that Tjrgg(Ti,±0 w ) = Oj so that the term involving this quantity may be 
dropped, from equation (36a). Substitution of equations (36) into 
equation (35) and. application of the previous expression (11) leads, 
for vanishingly small a, to the condition that 

f H t’ 0 (Tl,±e w ) dT) = 0 ( 37 ) 

-CO 

The boundary conditions ( 23 ), ( 24 ), ( 25 ), ( 34 ), and ( 37 ) are sufficient 
to determine the solution for ijr* in the holograph. 

As with f, the boundary-value problem for can be simplified 

from considerations of symmetry. Since if is antisymmetric with 
respect to 0, the nohhomogeneous boundary condition (23) which is 
imposed on if’ along the upper and lower boundaries must be symmetric 
in this variable. The remaining conditions, which are all homogeneous, 
are also symmetric. It follows that ty’ itself must be a symmetric 
function of 0 . 4 The problem can therefore be simplified by again 
el imina ting the lower half of the holograph and substituting in this 
case the condition that 

^©(tIjO) = 0 for !)■£- 1 ( 38 ) 

4 ' " ” ’ ’ "" 

This result can also be argued directly from considerations in the 
physical plane. It is necessary to make two observations as follows: 
(l) Since the profile itself Is symmetric about the chord line, the 
flow field at a negative angle of attack must be the inverted image 
of the flow field at an equal positive angle. ( 2 ) To be consistent 
with the basic perturbation assumption, it must be presumed that all 
changes in the flow field are smooth functions of angle of attack 
at a = 0 . These statements taken together imply that the vertical 
distance between any tiro points of equal q and corresponding posi- 
tive and negative 0 is, to a first order, unaffected by angle of 
attack. It follows that, for sufficiently small a, the increments 
in t and between the two points are equal and hence, on the basis 
of equation (19) , that’ the value of ty* at the two points is the same. 
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The problem which is finally to he solved is thus as summarized in 
sketch (e). The boundaries for this problem are identical with those 
used to obtain The boundary conditions are also identical insofar 

as the shock polar and sonic line are concerned. The only differences 
between the two problems are. in the conditions imposed along the bound- 
aries OB and OA. As was the case with the solution here must be a 
function of 0 W . 


0 



Because of the nature of the integral condition along the upper 
boundary OB, a direct solution for t' is not feasible by numerical 
methods. To obtain a solution, therefore, the problem is broken down 
into two subsidiary problems by means of the substitution 


t' = 'I'Y+^'b (39) 

where b is a constant whose value is to be determined. Boundary-value 
problems for ^’a and ijr'-g are then defined as shown in sketch (f ) . 
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In "both problems the integral condition along the upper boundary OB is 
ignored for the time being, and an arbitrary condition is introduced 
instead at the point E. In the problem for ^ , only the nonhomoge- 
neous condition ( 23 ) is imposed at the upper boundary, and the condition 
at E is the homogeneous one that =0* In the problem for f g, 

the homogeneous condition ty'g =0 is imposed along OB, and the condi- 
tion at E is that ty’-g has an arbitrary value £ 0* The conditions 

at the remaining boundaries are the same as in sketch (e) and are there- 
fore not repeated here. It is apparent that a superposition of 
and ty’jj will constitute a solution of the original problem provided the 
value of b is adjusted so that the integral condition (37) is satisfied 
on the upper boundary. The necessary equation for b is found by sub- 
stituting the expression (39) into condition (37) and is 


b 



/ 

-00 


O 

3 ♦W’b 0 *) dT l 


m 


Relations for quantities in physical plane. - To complete the funda- 
mental analysis, relations must be established between and the 

relevant quantities in the physical plane. Let x = x(t],0) and y = y(T],0) 
denote the coordinates at which a given velocity t ],0 is found in the 
physical plane when the profile is at zero angle of attack. As shown in 
reference 3 (pp. 29 - 31 ), the transformation equations (3), when applied 
to the case of zero angle of attack (and written in the present notation), 
can be put in the dimensionless form 

d f-f) = -i- (2T}f e dT) + V 0 ) ' (to) 


[(7+l) (t/c)] 1/3 d(j) = (^dT] + * 0 d0) = w (to) 

where Iw represents the integral 
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By talcing the origin of the physical coordinates at the leading edge and 
introducing the notation X = x/c and Y = [ (y+l) (t/c)] 1//3 (y/c ) , equa- 
tions (111) can he integrated to give 






,dri + t n de) 


Y(ri,e) = 


(20 v ) 1/3 


(^3a) 


(*3b) 


The integration in equation (43a) is performed in the hodograph over any 
curve C which begins at q = -00 and ends at the point q,0. The 
generalized coordinates X and Y at which the same velocity q,0 is to 
he found when the airfoil is at angle of attack are given correspondingly 

*>y 


X(r],e;a) = -k. P (2i]t e dT] + ijr de) 
«v Jq 

Y(r),0ja) = T|r 


(44a) 


(44b) 


The integration in equation (44a) is considered to he 'taken, over the 
same curve C as before. 5 The integral I w is now given by 


r° 

iw = / n^eCn^w-o-) 


It can he shown from equations (19) , (36), and (37) that for vanishingly 
small a 


^w _ 


Equations (44) can now he specialized in the light of the basic 
perturbation assumption. This assumption implies at once that the coor- 
dinates X and Y in the physical plane must he expressible in the form 


’if C lies slightly outside the domain in which is defined - as 
will he the case, for example, when the integration is talcen over the 
upper surface of the wedge in its undisplaced position - t{r is to he 
thought of as being continued analytically outside the boundary. 
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X(q,0;a) = X(q,0) + a X* (q,0) 

(47a) 

Y(q,0;a) = Y(q,0) + a Y* (q,0) 

(47b) 

where X* (q,0) =X a (q,0;O) and Y* (q,0) = Y a (q.0;O) . If expressions (19) 
and (47) are substituted into equations (44-), and equations (43) and (46) 
are taken into account, the following relations are finally obtained 
for X’ and Y* in terms of ijr': 

X ’ (T1,0) = 4 ^ f + % d0 ) 

(48a) 

(20 w ) 1/3 
r(n,e) =-~- — 

, «w 

(48b) 


The foregoing equations (48) give the initial rate of movement with 
angle of attack of a point of fixed velocity 11 , 0 . One requires for 
practical application, however, the rate of change of q and 0 at a 
point of fixed location X, Y. Equations relating the two sets of deriv- 
atives can he obtained as follows: If q and 0 are regarded in the 

physical plane as functions of X, Y, and ct - that is, q = q(X, Y; a) 
and 0 = 0(X, Y; a.) - then the corresponding total differentials are 

dq = q x dX + q-^dY + q a dct (49a) 

d0 > QjdX + 0ydY + 0 a da (49b) 

Consistent with the basic perturbation assumption, q and 0 can be 
written 

q(X,Y,a) = q(X,f) + a q‘(X,Y) (50a) 

0(X,Y,a) = 0(X,Y) + a 0> (X,Y) (50b) 

where q and 0 represent the conditions at a given point X, Y at zero 
angle of attack and q* and 0' are defined by q’(X, Y) = q a (X. Y; 0), 
0’(X, Y)=0 a (X, Yj 0). In view of equations (50), equations ( 49 ) can 
be written for vanishingly small a 

+ q^dY + q’da 
d0 = ©gdX + 0ydY + 0* da 



/ 


(51a) 

(511) 
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Similarly, from equations (4-7) , one can write for the differentials 
of X and Y as functions of T],0, and a 


dX = X^drj + X 0 d0 + X’da 


dY = Y^dn + Y e d0 + Y’da 

from which 


X^dTj + Xgd0 = dX - X'da (52a) 

YqdT) + Y 0 d0 = dY - Y'da (52b) 


Solution of equations (52) for dT] and d0 and comparison of the results 
with the alternative expressions (51) gives finally for and 0' 


, W' - * 0 Y ' 

V© - XgY^ 


(53a) 


•^T)^0 " XqYt, 


(53b) 


These equations can he put in more directly useful form by evaluating 
the derivatives of X and Y from equations (4-3) and substituting for 
X* and Y’ from equations (48). There results finally 


V(x,?) 


1 




(2Tit' e dTi + **40) 


- V’ 


(54a) 


0'(X,Y) ~—jL— \ J (2t)i|r^dq +*’,10) -2 t$ 0 ** 


(54-b) 


By means of these equations the initial rate of change of tj and 0 at 
some fixed point in the physical plane can be calculated corresponding 
to any chosen location in the holograph. The coordinates at which these 
derivatives apply are found from the solution at zero angle of attack 
by means of equations (4-3). 
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The foregoing equations are considerably simplified 'when applied at 
the surface of a wedge profile. Here the bound ar y condition is that 
is constant on a line of constant 0 (cf. equation (7))> with the result 
that if = 0. Equation (54a), for example, can thus be written as 
simply 

T) 

q‘(X,±0) i P dq (55) 

7te(n,±ew) J 


where the upper signs pertain to the upper surface and the lower signs 
to the lower surface. The corresponding rate of change of pressure 
coefficient is found by differentiating equation (4) with respect to 
angle of attack. If a is used now to denote the true angle of attack 
(related to the previously used, normalized angle of attack by an equa- 
tion like (lb)), such differentiation then gives 




(Va,-l) 1/a 


Here q’ is still the derivative with respect to the normalized angle 
as given by equations (54a) or (55)- With the aid of equations (5) and 
( 13 ), this result can be rewritten 


[(r+i)(t/c)] 1/3 (^) - - 2(20 w ) 1/s V 

atL y a=o 


(56) 


It can be seen from equation (55) and the symmetry properties of if 
and ijr* that q* must be of equal magnitude but opposite sign on the 
upper and lower surfaces of the profile. If the local lift coefficient 
is represented by Ap/qo= (Plover " Pupper) /lo* ft fb 611 follows from 
equation ( 56 ) that . 


[(r+l)(t/c)]Wa 


atAp/Q 

da 


= 4(20 w ) 1/3 q’ (X,+0) 


J a=o 


(57) 


where the notation q T (X,+0) indicates that the value is to be taken on 
the upper surface of the profile. Substitution from equation (55) gives 
finally 


[( 7 + 1 ) (t/c)] 


1/3 


jjAp/qp) 


L da 


-*a=o 


r\w,e,) ft, 

n* 0 (q,e w ) 


(58) 


By means of this equation, the initial rate_ of growth of lift at any 
chordwise station can be obtained. Since f and are both functions 
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of the parameter 0 W , the generalized quantity which appears on the left- 
hand side of equation (57) is also a function of this parameter. These 
results are in conformity with the rules for transonic similarity (see, 
for example, reference 11) . 


METHOD OF SOLUTION 


As in the previous calculations of ijr, the boundary- value problems 
for and fp can be solved through the use of finite -difference 

equations and relaxation techniques. A detailed description of the 
general method has been given in reference 3 and need not be repeated 
here. Most of the necessary finite -difference equations - notably the 
tedious ones along the shock polar and sonic line - can be taken over 
directly from the previous work. The only equations which need be 
altered are those directly influenced by the change in boundary condi- 
tions on the upper boundary and on the horizontal axis. The only real 
difficulty from this source is encountered in the solution for tJt’a in 
the vicinity of the shoulder (point B in sketch (f)). At the shoulder 
itself, the boundary conditions require a singularity in the first 
derivatives of t’A, which means that any purely numerical treatment 
would be of doubtful validity in the vicinity of this point. This dif- 
ficulty is overcome by subtracting out an analytical' solution of the 
proper singular form and then working locally with the difference 
between this solution and the desired unknown. The singular solution 
Is obtained from the general results of Guderley (inference 13) and is 
expressed in terms- of hypergeometric functions. The details of this 
and other matters regarding the numerical calculations for the front 
half of the profile are given in appendix A. 

With the solution known for the front half of the profile, the 
calculation of the lift on the rear half Is a simple matter. The com- 
putations are carried out in the physical plane and are based on the 
characteristics net previously constructed .for the flow over the rear 
wedge at zero angle of attack (see, for example, fig. 4 of reference 3) • 
Starting from the known solution for t’, one first employs equations (48) 
to compute the initial rate \ of movement of the points at which the Mach 
lines of the basic characteristics net meet the sonic line. Using these 
results and the known slope of the segments of the basic net, one then 
proceeds stepwise along consecutive downgoing Mach lines, calculating 
the initial rate of movement of successive intersection points on each 
line. By application of the proper boundary conditions at the surface 
of the wedge, the value of t] 1 at the surface is finally determined, and 
from this the initial distribution of lift is calculated. The details 
of the procedure are given in appendix B. 
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RESULTS AND DISCUSSION 


Calculations of the lift have "been carried out, following the 
methods just outlined, for the same values of Q v used in the work 
at zero lift, namely, 1.3, 1.6, 2.4, and 4.2. These values correspond, 
respectively (see equation ( 13 )), to values of the similarity param- 
eter i Q of 1,058, 0 . 921 , 0.703, ami 0.484. 

To illustrate the results for the front wedge in the hodograph, 
figures 1 to 3 have been prepared showing the variation of ijr’-g, 

and ty* for % = 1.6 . 6 The results for ijr'-g (fig. l) are only slightly 
different from those previously shown for f in figure 3 of reference 3 * 
As before, a rapid (hut regular) variation is apparent in the dependent 
variable in the vicinity of the point ^ = 0, 6 — 1. The results for 
(fig. 2) show *a rapid variation near the point q = 0, 6 = 0 W . This is 
a consequence of the previously mentioned singularity in the first 
derivatives of ijr’^ at that point. The values of i|f* (fig- 3) are 
found in the present case from the equation ijr* = - 0.5348 ijf'g 

(cf. equation (39))- They exhibit the same behaviob as does f in 
the vicinity of the singular point hut differ markedly in other parts of 
the field. For reference, the numerical values from which figures 1 and 2 
were plotted are given in tabular form at the end of the report. 

The complete results for the lift of the profile are given in 
figures 4 through 7- These results will he discussed in the following 
paragraphs . 


Chordwise Distribution of Lift 


Figure 4 is a plot of the calculated lift distribution, in tran- 
sonic similarity form, for the four values of | 0 considered in the 
present work. Also shown are the results for = 0 (Mq = l) given by 
Guderley and Yoshihara in reference 9- It is convenient for purposes 
of discussion to think of a similarity plot, such as that of figure 4, 
as pertaining to fixed values of t/c and 7 . From this standpoint, an 
increase from zero in the similarity parameter | Q is equivalent to an 
increase from 1 in the free -stream Mach number Mq. For simplicity, 
this point of view will sometimes be adopted in the descriptions that 
follow. 


®For the calculation of i|r’-g in this example, use was made of 236 

lattice points distributed as shown for ijr in figure 2 of reference 3* 
For t'*, 380 points were used with a distribution appropriate to the 
altered behavior of the dependent variable. 
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The lift distributions of figure 4 are al l , of the same general 
shape. In all cases the calculated lift 'tends toward infinity at the 
leading edge of the profile. This type of besult, which is of course 
physically impossible, is well known from the linear theory of airfoils 
at subsonic speeds. It is a result of the obvious failure of the small - 
disturbance approximations to conform with the actual phenomena in the 
vicinity of the leading edge. This local failure of the theory is 
known in the linear, subsonic case to be of little consequence insofar 
as the over-all lift is concerned. It may be presumed that a similar 
situation exists here. 

As one proceeds rearward from the leading edge, the lift distribu- 
tion falls more or less rapidly, reaching a value of zero directly 
forward of the shoulder. This latter result could have been foreseen, 
since the speed on both the upper and lower surfaces has a fixed 
(i.e., sonic) value at this location. Directly to jfche rear of the 
shoulder, the lift distribution starts anew from zetro. This must 
obviously be the case, since the expansion from sonic speed is, in 
Prandtl -Meyer flow, a unique function of the local turning angle, which 
is the same for both surfaces. Rearward from the shoulder the lift 
increases monotonically to a relatively small, finite value at the 
trailing edge. 

Over the front wedge, the four curves of the present study exhibit 
a uniform progression with respect to The curve of Guderley and 

Yoshihara, however, crosses the present curves at several points. The 
reasons for this are not clear, though it is highly unlikely that such 
a result could be in fact correct. The observed behavior may be due to 
some consistent inaccuracy in the present numerical approach or to the 
approximations introduced by Guderley and Yoshihara in satisfying the 
boundary conditions for the interdependent portion of the supersonic 
expansion fan. Over the rear wedge, the present computations give vir- 
tually a single curve for the four values of | Q . There is again, how- 
ever, a small inconsistency with the results given by Guderley and 
Yoshihara. This is as might be expected if the calculated flow over the 
front wedge is in error in either case. 


Lift -Curve Slope 


Figure 5 shows the generalized slope of the lift curve at zero 
angle of attack plotted as a function of the transonic similarity param- 
eter. Results obtained on the- basis of the transonic small -disturbance 
theory are shown by three solid-line curves. Each of these curves con- 
sists of two segments separated by a gap within which the curve cannot 
be defined on the basis of the available result's. The uppermost of the 
three curves gives the lift of the complete profile; the other two show 
the division of lift between the front and rear wedges. 
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The left-hand se gme nt of each of the curves in figure 5 shows the 
variation of lift-curve slope over most of the range of flight speed in 
which the how wave is detached, which is the range of pr im a r y concern 
in the present analysis. The calculated points from which these curves 
were drawn are shown in the figure. Hie points denoted hy squares were 
obtained hy mechanical integration of the lift-distribution curves of 
figure h . 7 The circled points on the vertical axis were located on the 
basis of the work of Guderley and Yoshihara. 

The right-hand segment of the curves in figure 5 shows the varia- 
tion of lift-curve slope in the range of flight speed in which the bow 
wave is attached and the flow is completely supersonic. To the order 
of accuracy of the present theory, this condition exists for the double- 
wedge profile at zero angle of attack when 2 1//3 = I.260. 8 Above 

this value, results completely consistent with the fundamental assump- 
tions of the transonic small-disturbance theory can easily be obtained 
by analytical methods. To this end, one need only presume that the 
speed is constant on each straight-line portion of the airfoil surface, 
a condition which is actually fulfilled over most of the pertinent 
range of f^. The necessary procedures are outlined in appendix C. To 
the accuracy of the transonic small -disturbance theory, the results 
provide an exact solution for the lift-curve slope of the front wedge 
for all values of £ 0 in the range of completely supersonic flow . For 
the rear wedge - and hence for the complete profile - the solution is 
exact down to a limiting value of | Q somewhat greater than 1.260. 

Below this limit the interaction of the shock wave from the bow and the 
expansion fan from the shoulder influences the flow over the rear wedge, 
with the result that the condition of constant speed is not satisfied. 
The position of this limit is difficult to determine exactly. As shown 
in appendix C, however, it must lie at a value of £ 0 less than I.287. 
The curves for the rear wedge and complete profile are thus approximate 
for at least a portion of the interval from 1.287 to 1.260 and" are 
therefore shown dotted in this range. It can be demonstrated that 
inclusion of the interaction effects in the analysis would cause an 
increase in the computed lift for the rear wedge. Exact results would 
thus lie somewhere above the dotted portion of the curves in figure 5. 


7 As in the earlier calculations of the drag coefficient at zero angle 
(cf. page 36 of reference 3), the integration over a small interval 
near the leading edge was carried out analytically on the basis of an 
asymptotic representation of the solution in the holograph plane. 

8 Attac hme nt of the wave takes place at the somewhat lower value of 

g 0 = 3 /( 1 j.) 2/3 = 1.191. For 1.191 < | 0 < 1.260 the wave is attached 
but the flow behind it is still subsonic. 
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The most interesting aspect of figure 5 is the behavior of the lift 
in the vicinity of shock attachment. Despite the gap in the curves in 
this vicinity, it is obvious that the lift-curve slope of the couple te 
profile must attain a maximum somewhere in the range from £ 0 = 1.056 
to £ 0 = I.287. This is in marked - and somewhat surprising - contrast 
to the previous results for the drag coefficient at zero angle of attack, 
which was found (reference 2) to decrease monotonically as the similarity 
parameter increased above zero. The peak in the curve in the present 
case is accompanied by a similar variation in the lift-curve slope of 
the front wedge. The results for the rear wedge may or may not pass 
through a minimum in the same range of i Q . 

A determination of the exact shape of the curves in the vicinity of 
shock attachment is not feasible on the basis of the present laborious 
methods. The existing curve for the complete profile does show a maximum 
in the range of completely supersonic flow, but this is in the portion 
of the range in which the computed curve is known to be erroneously low. 
If exact results were available for all values of £ , the maximum would 
undoubtedly be somewhat higher and displaced somewhat to the left. The 
infinity which appears in the slope of the curve at £ 0 = 1.260 (see 
appendix C) would probably disappear as well. The lift of the rear 
wedge, which now goes to zero at i Q = 1.260, would presumably remain 
finite throughout. 

Within the transonic range itself, the curves of figure 5 show 
little variation for some distance above a similarity parameter of zero. 
This is in accord with Guderley’s recent analytical study of two- 
dimensional flows with a free -stream Mach number close to 1 (refer- 
1 ence l4) . The results of Guderley’s work imply that, to the accuracy 
of the smal 1 -disturbance theory, the lif t-curve slope does not vary as 
the free -stream Mach number passes .through 1. The curves of figure 5 • 
have been faired so as to conform with this requirement. It is apparent 
that Guderley's result of zero variation may be taken as a good working 
approximation even at Mach numbers some distance removed from 1. The 
same result was found in reference 3 with regard to the drag coefficient 
of the complete profile at zero angle of attack. 

Over most of the range of completely supersonic flow, the lift- 
curve slope of the complete profile exhibits the type of variation well 
known from linear theory. This latter theory gives for the lift -curve 
slope of all thin profiles 


tel = k . 

(Mo 2 -l) 1/2 


(59) 
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which can "be ■written in terms of the transonic similarity variables as 


[(r + D(t/=)] 1/s ' (60) 

5 o 


The dashed curve in figure 5 is based on this equation. There is con- 
siderable quantitative difference between the 'linear and no nlin ear 
results for values of ^ just above 1-287- As' i 0 increases, how- 
ever, the curves given by the two theories appear to converge. This 
latter behavior is in accord with Spreiter's considerations (refer- 
ence 11 ) regarding the basic relationship between the linear and non- 
linear theories. 

To put the results in more familiar form, the lift-curve slope of 
the complete profile has been replotted in figure 6 as a function of 
Mach number for 7 = 1.4. The results of linear theory give a unique 
curve defined by equation (59)- The nonlinear, transonic theory pro- 
vides a family of curves with thickness ratio as a parameter . 9 As 
would be expected, the range of Mach numbers over which the linear 
theory is a poor approximation becomes smaller as the thickness ratio 
is reduced. It can be reasoned, in fact, that the nonlinear results 
must tend toward the results of the linear theory as t/c — > 0 . 


Center of Lift 


Figure 7 shows the chordwise position of the center of lift (x/c )i 
as a f unction of the transonic similarity parameter. The arrangement 
of the figure parallels that of figure 5* As before, the indicated 
points were calculated on the basis of the lift distributions of 
figure 4. The curve in the range of completely supersonic flow 
( 2:1.260) was obtained by means of the equations of appendix C. 

Only results for the complete profile are shown. 


g 

Here, as in reference 2, a multiplicity of curves could be obtained 
for each thickness ratio by using expressions for the pressure coef- 
ficient and similarity parameter different from those of equations (4) 
and (13) • In view of the recent developments outlined in footnote 2 , 
such complications now appear to be of lessened significance. 
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The movement of the center of lift with increasing Mach number is 
of some interest- At a free -stream Mach number of 1, the results of 
Guderley and Yoshihara indicate a position about 29 percent of the 
chord aft of the leading edge. As the Mach number is increased, the 
center of lift first moves forward, slowly in the initial stages and 
then more rapidly as the condition for shock attachment is approached. 
In the completely supersonic range, this trend is reversed; the center 
of lift then moves aft toward the midchord location given by linear 
theory. Apparently, the reversal of the direction of motion must take 
place rather suddenly in the vicinity of shock attachment. The limit 
of forward movement cannot be specified, except to say that it must lie 
somewhere ahead of 22 percent of the chord (and probably aft of the 
leading edge). The dotted (i.e., inexact ) portion of the curve passes 
precisely through the quarter-chord point at £ 0 = 1.260. (The corre- 
sponding lift distribution is one of uniform lift on the front wedge 
and zero lift on the rear.) Because of the interaction effects 
previously discussed, an exact result would lie somewhat above the 
dotted curve . 


CONCIHDMG REMARKS 


The present calculations add support to the growing conclusion 
(see references 2, 6, 7, 8, and lh) that no marked changes take place 
in characteristics of airfoil sections as the free-stream Mach number 
passes through 1. The establishment of this conclusion must be regarded, 
in fact, as one of the major successes of recent research in transonic 
flow. In the present case, as in the previous study of the drag coef- 
ficient at zero lift, the variation of the aerodynamic quantities with 
free-stream Mach number is most rapid in the vicinity of shock attach- 
ment. Unlike the behavior of the drag coefficient, however, the varia- 
tions here are large and characterized by a sudden reversal in the sign 
of the derivative. In drawing conclusions from these results it must 
be remembered, of course, that the theory assumes an inviscid medium 
and an airfoil of small thickness and infinite span. It also assumes, 
in effect, that the angle of attack is of an order smaller than the 
thickness ratio. To what extent the results will be valid for viscous 
flows about finite -span airfoils at practically usable values of the 
thickness ratio and angle of attack is difficult to say. The effects of 
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finite span, for example, will surely cause a reduction in the variations 
near shock attachment. In the present state of theoretical development, 
the study of these effects is a task for experiment. 


Ames Aeronautical Laboratory 

National Advisory Committee 
Moffett Field, Calif., 


for Aeronautics 
Aug. 1, 1952 
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APPENDIX A 

SOLUTION OF BOUNDARY-VALUE PROBLEM 
FOR FRONT WEDGE IN HODOGRAPH PLANE 


The solution of the boundary- value problems for and V'-g was 

accomplished by finite-difference methods similar to those developed for 
the calculation of V in reference 3* The description here will be 
limited to the few features wherein the present work departs from that 
discussed in the earlier paper'. (See general remarks under METHOD OF 
SOLUTION *) The notation and sketches follow the -conventions used in 
reference 3* 


Finite -Difference Equations Common to Both Problems 


The only finite -difference equations common to the problems for 
♦ ’A and V 'g but not found in the problem for ? derive from the bound- 
ary condition on the horizontal axis (see sketches (e) and (f)). This 
condition is given for both problems by equation (38) and is = 0 

for T] 5-1. In the previous work, the fin i te-difference equations for 
lattice points located on a boundary were obtained by approximation to 
the boundary condition itself. . In the present case, the approximation 
to the differential equation will be employed, and the boundary condition 
incorporated through use of the equivalent symmetry property. 

Consider a typical point 0 on the horizon- 
tal axis as shown in sketch (g). Point 3 is a 
fictitious lattice point located below the hori- 
zontal axis at 9 = -A, where A is the lattice 
interval. The finite-difference approximation to 
the differential equation (20) of the present 
text is given by equation (20 ) of reference 3 as 

v z + - 2T) 0 (r 1 + r 3 ) -2(1-2 ^ ' 0 = o (ad 

where tj 0 is the abscissa of point 0. The symme- 
try property leading to the boundary condition '(38) 

requires that t ' 3 = t’i, so that for points on the horizontal axis equa- 
tion (Al) reduces to 

r z + r 4 - Vi - 2 (i-2ti 0 )*' 0 = 0 


A Q\ 


I 

3 




(A 2 ) 
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The point at the intersection of the horizontal axis and the shock 
polar needs special consideration. Sketch (h) shows conditions at this 

point. Here, as before, point 4 is a ficti- 
tious point located below the boundary sym- 
metrical to point 3- It follows from the 
boundary conditions ( 25 ) and (38), both of 
which must be satisfied at the point 0, that 
the first derivatives in the coordinate 
directions are both zero at that point. On 
the basis of this fact, if the function 
ty'(q,0) is expanded in a two-dimensional 
Taylor's series about point 0, the following 
finite -difference relations for the second 
Sketch (h) derivatives are easily obtained: 

A 2 ir» I _ kilr> _ i -\fr» _ 1 ijfi 

^ “ qq|o l 2^ 2 2^° 

A 2 *' eelo = 2**3 - 2*' c - 

Here the symmetry property about the horizontal axis haB been used to 
equate to ^'g. Substitution of these relations into equation (20) 

for q = -1 leads to the following finite -difference equation for the 
point 0: 



4(1-2^)+^ - \ (l-2k 2 )*' 2 + 4* 


I _ 


4+1 (1-21C 2 ) 
2 


♦ ’ o = 0 


(A3) 


Finite -Difference Equations Special to 




B 


The only finite-difference equation special to the problem for ty’p 
is the one used to terminate the field of computation at some vertical 
line on the left. As in the corresponding work for this equation is 
derived from an asymptotic solution of the boundary- value problem valid 
for large negative values of q. The derivation is parallel to that 
described in detail on pages 16 and 17 of reference 3» 


%W=o 


'^■*0 
for eo 


ff g (y,0)=° 


0=0 W 



The boundary conditions which 
must be satisfied by ty'p at large 
negative values of q are shown in 
sketch (i). The shaded section show 
the anticipated variation of 
for constant q. A solution of the 
differential equation which satisfies 
the given boundary conditions is 


Sketch (i) 


0=0 
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+'8(1,0)= £ Cn cos(^)xy^K 1/3 [^ (-2,) 3 / 2 _ 
n=i 

■where K1/3 is the mod i fied Bessel function of the second kind of 
order 1/3 and the Cq are appropriate constants. If only the leading 
term of this solution is used and the Bessel function is replaced by the 
first term of its asymptotic expansion, there results 


P 


B (n,e) - c =os (J0 x(-i)- 1/4 ^ (-2i) 3/2 _ 


As in the earlier -work, let A denote the lattice interval and 0 some 
large negative value of q such that A/0 « 1 . It then follows from the 
foregoing solution that, to a first order and for a given value of 0, 


t’ B (-p-A,e) 



m 


By substituting this relation into equation (Al), a finite-difference 
equation can be obtained which is valid for points on the_ line q = -0 
and does not include any points to the left of this line (cf. equation (22) 
of reference 3 ) • 


Finite -Difference Equations Special to 


The only equations special to the problem for arise as a con- 

sequence of the condition along the upper boundary, where the values 
of are prescribed as a function of q. Along most of the boundary, 

this condition can be met by substituting the prescribed values directly 
into finite -difference equations of the type (Al) for’ points one interval 
below the boundary. Because of the nature of the boundary values 
near q = 0, however, some change from previous procedures is necessary 
in the vicinity of the shoulder. Modification is also required in the 
equations used to terminate the field on the left. 

Points near shoulder of wedge .- From the known behavior of f la 
the vicinity of the shoulder (see Guder ley's results, reference 13, for 
the flow around a convex corner), it can be shown that the variation 
of along the upper boundary near q = 0 must be of the form 

t' = *e(q,M -d H) 172 (A5) 

where D is a constant of proportionality. A singular solution of the 
differential equation (20) which is valid in the vicinity of the shoulder 
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and which satisfies the boundary condition (A5) is also obtainable from 
Guder ley's results. Tbis solution is, in the present notation. 


*' As (l.S) = I>(-’l) 1/2 (l-e) 1/S F^-i, i, i} jL) 


(A6) 


where F is the hypergeometric function and £ = £(q,0) is defined by 

, _ 9 ( 0 v -9) 2 

8 q3 

Equation (A6) is suitable for use near the upper boundary (0 = 0y> £ = 0). 
Near the sonic lin e (-tj=0, £ = - " ) the f ollowing alternate form is 
available: 





D 

2 2/3 



+ f (Qv-ef 


lVe 


F 


112. _1_\ 
6* 3* 3* 1 -t) 


(-n) 1/2 F fi 2 k. 1 M 

2 2 / 3 (i-^) 1 / 6 \6* 3' 3 J l-c/j 


(AT) 


If equation (A7) is evaluated on the sonic line, there results 



3 1 / 3 

27/e 


D (0 W -0) 1/3 


(A8) 


This result is in agreement with equation (32), which was developed from 
other considerations. It is apparent from equations (A5) and (A8) that 
a solution for ^'a will have a singularity in the first derivatives at 
the point q = 0, 0 = 0^. 

Because of the foregoing singularity, a direct numerical calculation 
of ijf'A might be expected to run into difficulty in the vicinity of the 
shoulder. Attempts along these lines lead, in fact, to the unlikely 
result of negative lift over a small region of the profile just forward 
of the midchord. Reductions of the lattice interval to quite small values 
served merely to decrease the extent of this region. This is in contrast 
to the previous work for f (and for as well), in which the singu- 

larity at the shoulder appears in the second derivatives. In that case, 
a sufficiently accurate solution for the unknown function could be obtained 
by direct calculation. In the present work, it was found necessary to 
subtract out the singularity in the first derivatives according to the 
following procedure: 

Let a function ST'a be defined such that 


5 *'A = t'A-V'Ag 
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where ty'^ is a singular solution of the type given by equations (a 6) 
and (A7). 8 If the actual, numerically determined values of ty'A on the 
upper boundary are examined, it is found that for a small length of the 
boundary near the shoulder these values can be replaced to a good approxi- 
mation by a 1/2 -power variation of- the form given by equation (A5) • This 
is done, and the constant D is determined such that within this length 
of boundary ty’« ^ ^w) or = 0. On this basis, a 

boundary- value problem for Sty ’a can be defined for a small region near 
the shoulder as shown in sketch ( j ) . The 
problem for Sty’A within this region is 
solved jointly with the problem for ty 'a 
in the remainder of the field. The two 
regions are fitted together by the use of 
overlapping lattices, much as is done in 
the case of a graded mesh (see reference 15 ) • 

The only difference is that equation (A9) I 

must now be utilized to make the transition A A 
between the two lattices at all their common Sketch ( j ) 

points. It is seen, from sketch (j) that 

conditions for oty 'a on both the upper boundary and sonic line are 
Identical with the cbrresponding conditions for ty. The finite -difference 
equations for the calculation of Sty 'a can therefore be taken over 
directly from the previous work. 

As nearly as one can judge from experience with various lattice 
spacings, results obtained by the foregoing process are quantitatively as 
well as qualitatively reliable. The primary source of error is in replac- 
ing the actual values of ty T A along the upper boundary by a l/2-power 
variation. Since the region over which this Is done in the holograph 
corresponds to a very small portion of the chord in the physical plane, 
errors from this source are probably small . 



Points far to the left .- The 
boundary conditions for ty'A at 
large negative values of q are 
shown in sketch (k). From the 
asymptotic solution for the basic 
problem (equation (21) of refer- 
ence 3), the expression for ty'A 
along the upper boundary is found 
to be 



Sketch (k) 
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t’ A (Tl,ew) = 1 f e (n^) = E (-Ti) -1/4 exp 


- 2 L ( - 2 t |) 3 ^ 2 

38v 


(A10) 


where E is a constant. 

Because of the nature of the houndary condition (AlO) , it is not 
possible to write an asymptotic solution for ty’ A for large negative 
in a single term. For this reason, the procedure previously used to 
terminate the field of calculation at some location on the left cannot 
be applied in the present case. An alternative procedure, somewhat more 
arbitrary in nature, can be devised by writing * A in the form 


r A 0i,e) = t' A (Ti,e w ) + 5t’ A (T],0) 


(All) 


where B^’a is now defined by S^’a (tj, * A (q,0) ~ ^ 'a ( t 1,®w) ( see 
sketch) . The attenuation of ^’ A in going from a point at q = -p 
to a point at rj = -3 - A is then found by treating each of the terms 
in equation (All) as an independent quantity. The attenuation of 
y a ( t 1,0v) ts found from equation (AlO) by a procedure similar to that used 
in obtaining equation (At-) . The result is 


V A (- P-A,e w ) 






(A12) 


To obtain a corresponding equation for 8^' A , it is assumed that for a 
given value of 0 this quantity attenuates in the same manner as was 
previously found for ty’-g. One thus has from equation (A^) 


8*' a (-P-A,0) 




S*’ a (-P,0) 


(A13) 


Substitution of these expressions in equation (All) for 1 ] = -P 
gives finally 




= (i-$) 


)x 


| exp ♦’ A (-^^ + 


A 


(Al4) 


Since t’A^j^w) is a known quantity for any given value of p, this 
equation can be used to terminate the field of calculatidn in the same 
manner as was done with equation (Ah) . The considerable element of 
arbitrariness in the derivation of equation (Al4) can be tolerated since 
the over-all solution is insensitive to changes in the left-hand portion 
of the field. 
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Solution of Finite-Difference Equations 


The techniques used to obtain a solution of the finite -difference , 
equations for i}^ and ty'-g were the same as those described in refer- 
ence 3 for the basic solution of ijr. In general., the graded lattice as 
used for f (see fig. 2 of reference 3) was suitable for the' solution 
of For however, different gradations were necessary with 

the smallest lattice spacing being used near the shoulder (point B of 
sketch (e)). The value of t’g at the intersection of the shock polar 
and the sonic line was chosen as 10,000 so that the previously obtained 
values of ijr could be used to provide the . initial guess for tJt'b • 

In the course of the present work, a useful technique was found for 
locating regions of relatively large error in the numerical solution. 

By use of one form of Green's theorem plus the differential equation (20), 
it can be shown that around any contour enclosing a region in which 
equation (20) is satisfied the following relation must hold: 


/ 


(2Tl*»od7} + 


de) = 0 


(A15) 


In a numerical solution the line integral in equation (A15) will not, 
except by rare coincidence, be precisely zero around any given contour. 
The amount by which it differs from zero may be taken as a rough measure 
of the adequacy of the numerical solution over the region within the 
contour. If the entire field of calculation is subdivided into a number 
of contiguous regions, it is thus possible, by evaluating the integral 
around each of the enclosing contours, to locate regions within which 
the error is relatively high. The solution in these regions can then be 
improved by advancing locally to a finer mesh. This technique was found 
to be of great "help in the present work. It would probably be useful in 
other elliptic boundary-value problems for which a relation analogous to 
equation (A15) can be obtained. 



42 


NACA IN 2832 


APPENDIX B 

CALCULATION OF FLOW OVER REAR WEDGE 
IN PHYSICAL PLANE 


The procedure used to calculate the flow over the rear wedge has 
"been outlined in the section METHOD OF SOLUTION. The fundamental opera- 
tion is to determine^ hy stepwise methods, the initial rate of movement 
of the known intersection points in the basic characteristics net. The 
methods which are used depend on the fact that these points are, hy 
virtue of the basic characteristics construction, points of fixed T),0 
(cf . equations (55) and (57) of reference 3) • 

The first step is to determine the ini tial rate of movement of 
those points at which the Mach lines of the basic characteristics net 
meet on the sonic line. For this purpose, consider equations (48) , 
which give the initial rate of movement of a general point of fixed T},0. 
If these equations are specialized to apply to points on the sonic line, 
the following relations are obtained: 

x'(o,e) = 


Y'(o,0) = 

To write equation (Bla) the path 
taken along the upper boundary from 0 to B (see sketch (e) on page 21) 
and thence downward along the sonic line. The contribution of the 
portion from 0 to B is zero by virtue of the condition (37)- In apply- 
ing these equations, the value of is known from the basic solution. 

The integral in equation (Bla) is evaluated by mechanical integration 
of a curve of numerically determined derivatives. Proper allowance is 
made for the singularity at the shoulder by integrating the singular 
solution analytically. The component rates of movement of the sonic 
point at the shoulder are both seen to be zero. 

The next step in the solution is to calculate the rate of movement 
of intersection points downstream of the sonic line. This is done by 
proceeding stepwise along consecutive downgoing characteristics. 




d0 


(Bla) 


_(20j 

hi. 


1/3 




(Bib) 


w 


of integration in equation (48a) is 
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Consider three typical net points as shown in sketch (l) (cf. also 
fig. ^ of reference 3) • The dashed lines represent the original posi- 
tion of the Mach lines through points a. b, and c, and the solid lines 


represent their displaced posi- 
tions corresponding to a small, 
finite angle of attack a. 

Since the intersection points in 
the Mach net are points of fixed 
t the components of their 
displacement are given by aX 1 
and aY* . The slope of each segment 
of Mach line is taken, in accord 
with the procedures of refer- 
ence 3, as the average of the 
slopes calculated at the two end 
points. The slope calculated at 
each end point depends, fn turn, 
only on the value of tj at that 
point (cf. equation (5k) of 
reference 3) * 



It is desired now to deter- 
mine X f and Y* at point c in Sketch (l) 

terms of X 1 and Y* at points a 

and b. Since the value of T) at a given net point is the same in the 
displaced and undisplaced positions, it follows from what has been said 
above that each segment of Mach line must retain its orig in al slope 
after displacement. If this slope is denoted by m, the following 
relations are then readily obtained: 


X’c 


Y , a-Y*b+ttb c X t b-aacX , a 
“be “ ®ac 


(B2a) 


Y’ 


Pbc^'a "^acT'b+^ac^bcCX’b -X'a) 
®bc “ 


(B2b) 


With these relations, it Is a simple matter to calculate the initial 
rates of movement of successive net points on consecutive downgoing 
characteristics. For the first characteristic to be considered, point 
b is taken at the shoulder of the profile, where X 1 and Y r are both 

zero. Thus, X' c and Y’ c for net points on this characteristic can be 

determined solely in terms of X’ a and Y’ a and the slopes ma, c and £tb c . 

For the remainder of the downgoing characteristics, X’-fc and Y't, are 
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known from calcinations along the characteristic immediately preceding. 
The actual calculations can "be carried out in straightforward tabular 
form. 


The foregoing procedure enables the calculation of X’ and Y’ for 
all net points except the ones originally at the surface of the rear 
wedge. For these points, consideration must be given to the required 
boundary condition at the surface. This boundary condition is 

0(X,+O;a) = - (0-^+a,) (B3) 

from which it follows that 

0*(X,+O) = - 1 (b4) 

The problem now is to determine X’ and Y' at the surface of the 
wedge in such a way that equation (b 4) is satisfied. To do this equa- 
tion (53b) is first specialized to the surface of the wedge, where it 
is readily shown that Xq = Y^ = 0. In view of condition (b 4), there 
results 


Y’(T),-ew) = Ye^-Sw) 


(B5) 


The value of Y* at points originally on the surface of the wedge is 
thus fixed directly by the basic solution. The corresponding value of 
X’ can be found from a construction analogous to that of sketch (l) 
and is 


_ Y ’a-Y’ c -fiacX’a 

A c 

-mac 


(B6) 


The point c is now the point originally on the surface of the wedge 
(i.e., Y’ c is as given by equation (B5) ) , and the remaining notation 
is the same as in sketch (l) . 


Application of equation (B5) requires the knowledge of Yg(q, -0^ , 
which in the case of the wedge profile is equal to l/0y. Evaluation 
of the latter derivative can.be carried out directly from the basic 
Mach net, but the procedures are cumbersome and inaccurate. A better 
method is to use the equations of motion (cf. equation (6) of refer- 
ence 4) to express 0y in terms of Following this procedure, one 

obtains finally 




' (23w) 1/3 

25i(X,0)ii x (X,40) 


(B7) 
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The quantities fj and which appear here are easily evaluated from 
the basic solution for the chordwise distribution of rj. 

The preceding equations enable the calculation of the initial rate 
of movement X T for points originally on the surface of the rear wedge . 
The final step is to determine the corresponding distribution of lift. 
For this purpose, equation (53a) is specialized to points on the rear 
wedge to obtain 


V(x,40) = - x‘(T!,-e w )/x-(fi,-e w ) 


which, in view of the boundary conditions, can be shown to be equivalent 
to 


T)*(X,+0) = - X»(Ti,-ev)t} x (X,+0) (B8) 


The distribution of lift is then obtained 1 from equation ( 57 )- 
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APPENDIX C 

SOLUTION OF PROBLEM FOR COMPLETELY SUPERSONIC FLOW 


Calculation of Lift-Curve Slope arid 
Center of Lift 


If conditions are such that £ 0 5f2 1//3 = 1.260 (corresponding 
to c f. equation (13)), then the basic flow over the profile at 

zero angle of attack is completely supersonic. The solution for the 
lift-curve slope and center of lift at a vanishingly sma ll , angle of 
attack can then he carried out analytically as follows: 

Consider a completely supersonic flow about the double -wedge profile 
at a small angle of attack. In the physical plane the flow field has 
the we 11 -known appearance shown on the left in sketch (m) . The corre- 
sponding hodograph of the flow along the upper surface, in terms of the 

0 



normalized small-disturbance variables T) and 0, is shown on the right. 
The quantities 6 V and a are, as before, the half -angle of the wedge 
and the angle of attack (also normalized; . Except for a small range 
of i Q just above 1.260 (see below), flow conditions must be constant 
along each of the segments 1 and 2 of the upper surface. In the hodo- 
graph each of these segments is thus represented by a single point 
located as shown. It is apparent that for a given value of 0 W , the 
speeds and tj 2 , which are the primary unknowns in the problem, are 

functions solely o£ .the angle of attack a. 
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To find the lift-curve slope and center of lift it is necessary- 
first to find the derivatives T) ^ = (dT] 1 /da) a _ Q and rj =(dT] 2 /da) a _ 0 . 
This can he done with the aid of the equations for the shock polar 


Q = (l-t) / 1+Ti 


and for the downgoing characteristic 


8 = constant - — 1— t\ 3/s (C2) 

3 

To find T] 1 !, one must utilize the houndary condition 6 1 = 9^ - a. 
Substitution of this condition into equation (Cl) provides the following 
implicit equation for rj 


8^ - a = (1-Tjj) */l+'n 1 


Differentiation of this equation gives 

da 1+ 3 t 1 1 

dT 1l 2/1+71! 

From this it follows that 


V, = 


2/l+fj 1 

1+3^! 


where, as in the text, the bars denote the value of at a = 0. 

The value of ij can he found in terms of the parameter Oy by solving 
equation (C3) for iq x with a set equal to zero. The result, obtained 
through standard methods for the solution of cubic equations , is 


tu = - 


^3/2 9y 
2 cos 

3 


where 


<p = arc cos 


Compare equation (53) of reference 3« 
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To find t)* 2 , equation (C2) for the downgoing characteristic is 
first specialized so as to pass through the point 1. This gives 


e - (e v -a) + (^3/2 . ,,3/2) 

Substitution of the bou ndar y condition 9g = - 9y - ct. then provides 
the result that 


= (tl 


^2 _ \4i 


3/2 


\ 2/3 




J 


Taking the 'derivative with respect to a, one obtains finally at a = 0 


1/2 


Vo - 


\ 1/3 


Vi 


(c 6) 


( v/a+ vH . 

where is given by equation (C4-) and fj ■ by equation (C5) . 

Since the value of i] 1 is constant on each segment of the profile, 
the lift-curve slope is easily found from equation (57) and is 

f cLc \ 

[(7+l)(t/c)] V3 \-~\ = 2(20 V ) 1/ ' 3 (T 1 * 1+ T ] « a ) 

v 'a=o 


Substitution from equation (c6) gives 

[(7+l)(t/c)] V3 C^) = 2(2e w ) 1/3 i ) ' 1 

v '0=0 


« x / 2 


1 + 


(>.•'■> p -r 


(01) 


The moment-curve slope, for moments taken about the leading edge, is 
found to be 

[(r+i)(t/o)] 1/3 (~) - - | 3n' e ) 


or 

[(7+1) (t/ c)J 


1/3 ( ^ c m \ 


= - i (20 w -) 1/3 T)’ i 


1 + (* 


3^ 




/ * 


(08) 
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The position of the center of lift is given accordingly by 


± 3ij 1 1/g 

^ _ (dcm/aa)^ ! ^ (3//5)<Hr] 1/S 

' " (d^j/to)^ = 1 » X x + ij, 1 /* 

I T5i S/2 + (3-/2)8„] 1/3 , 


In equations (C7) and (c8), the first term inside the brackets repre- 
sents the contribution of the front wedge, the second term that of the 
rear. 


Equations (C7) and (C9) are the basis for the curves shown in 
figures 5 and 7 for values of £ o ^1.260. The results show certain 
curious features when the flow over the front wedge is just sonic 
(i.e., 7^=0, 0^=1, § o =l-260). These are as follows: 

(a) The lift contributed by the rear wedge is zero (see equa- 
tion (C7). 

(b) The center of lift is at the quarter-chord point (follows 
from statement (a) .plus the condition of uniform lift on the front 
wedge; see also equation (C9)). 

(c) The rate of change with respect to g is infinite both for 
the lift-curve slope of the complete profile ana for the position of 
the center of lift (follows from differentiation of equations (C7) 
and (C9) ) - 

These results are associated in every case with the behavior of the lift 
calculated for the rear wedge. 


Estimation of Lower Limit for Constant Speed 
Along Rear Wedge 


The features just enumerated, though having a certain curiosity in 
themselves, cannot be accepted as completely correct. Because of inter- 
action effects between the shock wave from the bow and the expansion 
fan from the shoulder, the fundamental condition of constant speed at 
the surface of the profile will not be satisfied along the rear wedge 
until the value of g Q is somewhat greater than 1.260. Until then, 
disturbances reflected from the shock wave will reach the rear wedge 
and cause a slight decrease in speed toward the trailing edge. This 
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effect will cease when the forwardmost reflected Mach wave just touches 
the trailing edge. The exact value of £ 0 at which this condition will 
he met is difficult to determine. An upper hound can, however, he esti- 
mated as follows: 

Consider the basic flow field (o=0) over the upper half of the 
profile when the first reflected Mach wave just strikes the trailing 
edge. Sketch (n) shows such a flow field as it would appear in tran- 
sonic similarity form (cf. pp. 12-13 of reference 2). In drawing the 
sketch a special assumption has heen introduced beyond those implicit 
in the small -disturbance theory; namely, that the first reflected Mach 



Sketch (n) 

wave is straight and has an angle of inclination p equal to that of 
the first wave in the expansion fan. With this assumption, the corre- 
sponding value of § 0 is easily determined. Since the reflected wave 
must actually be curved downstream, the value so determined will be 
greater than the correct value for the required condition. 

On the basis of sketch (n), the following equation can be written 
between the shock angle X and the Mach angle p: 

tan X = i tan p (CIO) 

3 

A relation between the shock angle X and the speed in the region 

behind the shock can be obtained from equation (C3) and the known prop- 
erties of the shock polar. The result is 
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tan X = 


= (29y) 

Vl+% 


1/3 


An analogous expression for the Mach angle 4 is given by equation (5*0 
of reference 3 and is 


tan (i = 


( 26 *) 


1/3 


</ 2 ^~ 


Substitution of these relations into equation (CIO) and solution for 
gives 

^■IF 


The accompanying value of Qy, found from equation (C3) with a=0, is 


0 v = O .9685 


This corresponds, according to equation (13), to 

£ 0 = 1.287 (Cll) 

Thus, for values of g between 1.260 and some limit less than 1 . 287 , 
the results of equations (C7) , (C8), and (C9) are not exact insofar as 
the contribution of the rear wedge is concerned. It can be reasoned 
that In this range an exact solution would indicate more lift for the 
rear wedge than does the present analysis. 


( 
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TABLE XI.- VALDES OF P A FOR. % = 1.6 (£ o =0.92l) 
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Figure 5. - Generalized lift-curve slope as a function of transonic similarity parameter. 
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